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Abstract. We analyze optical (UV rest-frame) spectra of X-ray selected narrow-line QSOs at redshift 1.5 ^ z ^ 3.7 
found in the Chandra Deep Field South and of narrow-line radio galaxies at redshift 1.2 ^ z ^ 3.8 to investigate the 
gas metallicity of the narrow-line regions and their evolution in this redshift range. Such spectra are also compared 
with UV spectra of local Seyfert 2 galaxies. The observational data are inconsistent with the predictions of shock 
models, suggesting that the narrow-line regions are mainly photoionized. The photoionization models with dust 
grains predict line flux ratios which are also in disagreement with most of the observed values, suggesting that 
the high-ionization part of the narrow-line regions (which is sampled by the available spectra) is dust-free. The 
photoionization dust-free models provide two possible scenarios which are consistent with the observed data: 
low-density gas clouds (nn ^ 10^ cm~^) with a sub-solar metallicity (0.2 <, Zgas/^© ^ 1.0), or high-density gas 
clouds (riH ~ 10® cm~^) with a wide range of gas metallicity (0.2 <, Zg^s/ZQ <, 5.0). Regardless of the specific 
interpretation, the observational data do not show any evidence for a significant evolution of the gas metallicity 
in the narrow-line regions within the redshift range 1.2 <i z ^ 3.8. Instead, we find a trend for more luminous 
active galactic nuclei to have more metal-rich gas clouds (luminosity-metallicity relation), which is in agreement 
with the same finding in the studies of the broad-line regions. The lack of evolution for the gas metallicity of the 
narrow-line regions implies that the major epoch of star formation in the host galaxies of these active galactic 
nuclei is at z 4. 
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• I— I 1. Introduction only for relatively bright targets (e.g., Teplitz et al. 2000; 

^ , , , . , , . . r Pettini et al. 2001; Teczaet al. 2004; Shapley et al. 2004). 

L_i Understanding galaxy formation and evolution is one of t ^ j i- i x- i ■ /a/-(at \ i -u-^ u • i i 

i2 , ,1 instead, active galactic nuclei (AGJNs) exhibit bright emis- 

^ , the key astropnysical issues which is being pursued in - r . , r rr^r i j-u t,-t, u j 

1, ., .. r i • sion lines at rest-irame UV wavelengths, which can be used 

this decade. The chemical composition of gas and stars in , . . ,^ ^ n- • i - i t„ • ^ 

, .... to investigate the gas metallicity even m mgh-2; objects, 

galaxies provides important information on this issue be- o ^ • 1^ i- r 1 • i r\c<r\ u- -L^^ j. 

° . , / . 1 . ... Spectroscopic observations of high-z QSOs suggest that 

cause it IS a tracer 01 the star formation history m galaxies. , n- -j. • ;i u j r ■ /t^t n\ j. j j. 

° gas metallicity m the broad-line region (I3LR) tends to 

A promising way to study the chemical evolution of galax- ^^^^ ^.^^^^ ^^^^ ^^j^^, U^mann & Ferland 1992; 

ies is to measure their metallicity as a function of redshift. t->- ^ • i ^ 1 onno at ^ i nnnr\ i- 1 

/ . Dietrich et al. 2003; JNagao et al. 2005) reaching as much 

Since it is extremely hard and time-consuming to measure ^ 1 r ^ 1 j • ^ 1 nnnn\ tt ■ i 

„ „. . r r ■ 1-1 1 ■ 1 1-1 as Z ^ 15^5© (Baldwin et al. 2003). However it is not 

stellar metallicity of taint high-z galaxies, because high , , j. n- • r j r j r 

1 clear how the gas metallicity interred from the broad lines 

quality spectra 01 shallow absorption features are required . -l - i r^i^u^ i- 

7, . • 1 r-ri^^Ti,^,, 1 r^ri^^N • IS rclatcd to thc chcmical propcrtics ot thc host galaxics, 

(but see, e.g., Pettim et al. 2000; Mehlert et al. 2002), m- . ,1 1 j r r a/-iat 1 i n 

^ , 1 , , since the broad lines of AGNs sample only a very small 

vestigatmg the gas metallicity through emission lines IS a . r i 1 • /n ^^1 • j. 1 

. . °. , region of galactic nuclei (-Kblr ^ 1 pc; e.g., Kaspi et al. 

promising strategy. However, most of the available diag- r,nnn\ u- 1, i, 1 j -ji -ti, 4-u 1, 4. 

. . . . ' . , . „ 2000), which may have evolved more rapidly than the host 

nostic emission lines associated with massive star forma- 1 

tion (e.g., [Oii]A3727, [Oiii]AA4959,5007, [Nii]A6548,6583) , 

, r , • 1 1 ,1 1 ;i 1 J An alternative possibility is to use narrow line AGNs. 

are m rest-trame optical wavelength and thus are shitted ^ ,.a^ 

., ■ r 1- A J-1U In this paper we focus on narrow- line AGNs at high red- 

mto near mtrared m high-z galaxies. Accordingly, observa- in 

,. r;i .... • 1,- u 1 • f -ui shift, in particular on high-z radio galaxies (HzRGs) and 
tions 01 these emission hues m high-z galaxies are feasible ' JL^ x S., , . . .^^^ 
^ ^ type-2 QSOs (QS02s). The nuclei of narrow-Une AGNs 

Send offprint requests to: T. Nagao are thought to be obscured by edge-on optically thick tori 
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(e.g., Antonucci & Miller 1985; Antonucci 1993; Cohen et 
al. 1999). Since the broad emission lines and the strong 

ionizing continuum arc blocked by this "natural coron- 
agraph" , we can investigate narrow UV emission lines 
whose spatial extension is roughly comparable to that of 
the host galaxies 10^^* pc). without any complex de- 
blending of broad and narrow components for the emis- 
sion lines. Although narrow emission lines of HzRGs are 
often strongly influenced by a radio jet in terms of kine- 
matics and morphological properties (e.g., McCarthy et 
al. 1991; Baum & McCarthy 2000), the ionization mech- 
anism is mostly dominated by photoionization, not by 
shock ionization (e.g., Villar-Martin et al. 1997; Allen et 
al. 1998; Iwamuro et al. 2003). Therefore, we can obtain in- 
formation on the gas metallicity of the narrow- line regions 
(NLRs) by comparing the observed emission-line flux ra- 
tios with the predictions of photoionization models. By 
comparing the NLRs of HzRGs and QS02s with low-z 
type 2 AGNs (Seyfert 2 galaxies; Sy2s), we can investi- 
gate whether the gas metallicity evolves on a cosmological 
timescale or not. 

By focusing on the flux ratio of NvA1240/CivA1549 
that is one of the most frequently used metallicity diag- 
nostics for AGNs (generally for the BLR; e.g., Hamann 
& Ferland 1992, 1999; Dietrich et al. 2003), De Breuck et 
al. (2000) reported gas metallicity of HzRGs in the range 
OAZq <^ Zgas ^ 3.0^0 or possibly even much higher (see 
also van Ojik et al. 1994; Vernet et al. 2001). They also 
claimed a metallicity evolution within their sample from 
z > 3 to z < 3: the gas metallicity of all HzRGs at z > 3 in 
their sample is Zgas < '2Zq, at variance with their sample 
at z < 3. Norman et al. (2002) reported that the strong 
NvA1240 emission of a X-ray selected high-z QS02 fomid 
in Chandra Deep Field South (CDFS; Giacconi et al. 2002; 
Rosati et al. 2002), CDFS-202 (z ~ 3.7), is consistent with 
a super-solar metallicity of its NLR, and more specifically 
Zgas ~ 2.5 — 3.0Zq. However, since the NvA1240 emission 
becomes very weak for metal-poor gas (i.e., Z^^s ^ ^Zq), 
only upper- limit fluxes on NvA1240 are available for the 
majority of the HzRG sample in De Breuck et al. (2000), 
which makes the investigation of the metallicity evolution 
difficult. An additional issue is that Iwamuro et al. (2003), 
based on rest-frame optical spectra obtained by sensitive 
near-infrared spectroscopy, recently reported that HzRGs 
at 2.0 < z < 2.6 tend to show sub-solar metallicities 



gas 



O.2Z0), significantly lower than values reported 
by De Breuck et al. (2000). Although it is clear that im- 
proved and additional observational data for a larger sam- 
ple of HzRGs are required to reconcile this disagreement, 
both more sensitive measurements of NvA1240 and deeper 
near-infrared spectroscopic data are very tough to obtain 
for HzRGs and QS02s. 

To tackle these issues, it is useful to find alternative 
metallicity diagnostics that consists only of strong UV 
emission lines. In this paper, we investigate strong UV 
emission lines for HzRGs, QS02s, and Sy2s. By compar- 
ing the compiled data with photoionization model calcu- 
lations, we discuss the evolution of the gas metallicity of 
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Fig. 1. Redshift distributions of the sample objects, (a) 
The HzRG sample. The white histogram denotes the orig- 
inal HzRG sample (at 1.2 < z < 4.0; 114 objects) of De 
Breuck et al. (2000) while the black histogram denotes 
the 49 HzRGs whose emission-line fluxes of CivA1549, 
HeiiA1640 and Ciii]A1909 are measured, (b) The QS02 
sample. The white histogram denotes the AGNs (includ- 
ing both broad- line and narrow- line AGNs) at z > 1.2 
found in the CDFS (51 objects; Szokoly et al. 2004) and 
one additional QS02 reported by Stern et al. (2002), while 
the black histogram denotes the narrow-line AGNs among 
them (10 objects) which are presented in Table 1. 



the NLRs in AGNs, from high-z to the local universe, us- 
ing only strong UV emission lines. The compiled data and 
our photoionization model calculations are presented in §2 
and §3, respectively. We compare the observational data 
with the model results, discuss the properties of NLR gas 
clouds and the implications for the chemical evolution of 
AGNs in §4. A brief summary is given in §5. 
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Table 1. Emission-Line Measurements of X-ray selected QS02s at 2; > 1.2 
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Object Redshift LyaA1216^ NVA1240'' CivA1549^ HeiiA1640^ Cin]A1909'' 



CDFS-027 


3.064 


12.6±0.7 


2.5±0.7 


6.4±0.5 


2.3±0.9 


<2.9'^ 


CDFS-031 


1.603 






24.1±1.4 


13.3±1.2 


10.3±1.3 


UJJr b— Uo / 


Z.ODZ 


110 0_L1 Q 


o.4±1.4 


1 / .o±U.8 


( .D±U.O 


io.o±U.9 


CDFS-112a'' 


2.940 


175±6 


22±6 


30±6 


19±6 


15±6 


CDFS-153 


1.536 






25.5±1.4 


6.2±1.1 


13.7±1.6 


CDFS-202 


3.700 


78.1±1.0 


26.8±1.1 


38.9±1.2 


19.7±1.5 


<12.9'* 


CDFS-263b 


3.660 


70.9±0.7 


4.6±0.7 


15.5±0.8 


<4.0'^ 


<7.6'* 


CDFS-531 


1.544 






22.0±1.4 


17.4±1.5 


14.4±1.5 


CDFS-901 


2.578 


37.1±0.6 


6.5±0.8 


19.7±1.0 


<2.8'* 


3.3±0.9 


CXO 52" 


3.288 


189±4 


6±1.2'' 


35±2 


17±2 


21±2 



^ Fluxes are in units of 10~^* ergs s~^ cm"^. 

FITS file not available for this source at the CDFS web site at the time of writing of this paper; the fluxes and their errors 

were estimated from GIF file on the web site. 
" Data taken from Stern et al. (2002). 

3cr upper-limit flux. 

The error was assigned by us; see text. 



2. Data 

2.1. Data compilation 

To investigate the possible chemical evolution of the NLRs 
in AGNs, we compiled the fluxes of some strong UV emis- 
sion lines, NVA1240, CivA1549, HeiiA1640 and Ciii]A1909. 
We focus only on these emission linos since fainter lines 
were measured only in a few type 2 AGNs. The details 
of the data compilation for each population of objects are 
given below. 

(1) HzRGs: We used emission-line fluxes of HzRGs 
given by De Breuck et al. (2000), who compiled emission- 
line flux data of 165 radio galaxies from the literature. This 
sample contains 114 HzRGs at 1.2 < z < 3.8 for which we 
can investigate the rest-frame UV spectra. We assigned 
the error on each line flux by checking the original refer- 
ences. For objects whose emission-line flux errors are not 
given in the original reference, we assigned a 20% error for 
each emission line. For emission-line fluxes measured with 
a statistical significance lower than 3cr, we adopted 3cr 
upper-limit values instead of the measured fluxes. We also 
referred to the measurements of omission-line flux data 
of 9 HzRGs presented by Vernet et al. (2001). Although 
these HzRGs are included in the sample of De Breuck et 
al. (2000), we adopted the data of Vernet et al. (2001) 
of these 9 HzRGs because of the higher signal-to-noise 
ratios of the Vernet et al. (2001) data. Accordingly, the 
number of HzRGs whose fluxes of CivAl.549, HeiiA1640 
and Ciii]A1909 are measured is 51 in total, of which two 
objects are lower- or upper-limit data. 

(2) QS02s: We used the spectra of X-ray selected 
QS02s in CDFS recently obtained by Szokoly et al. 
(2004). Among 288 objects whose spectra were obtained 
with VLT/FORS by Szokoly et al. (2004), 51 objects are at 
z > 1.2. From these 51 objects, we selected 9 objects with 
a high quality spectrum that shows only narrow emission 



lines. All of these 9 objects (distributed at 1.54 ^ z ^ 3.70; 
see Table 1) show strong CivA1549 emission which in- 
dicates, along with their X-ray emission, that these ob- 
jects harbor an AGN (i.e., not starburst galaxies). We ob- 
tained spectral data of these 9 objects from the CDFS web 
site^, and measured their emission-line fluxes and errors by 
means of a simple Gaussian fitting. The measured fiuxes 
are given in Table 1. In this table fluxes of Lya are also 
given for the reader's convenience, although we do not use 
the Lya flux in the analysis and discussion in this paper. 
The reported errors do not take any possible systematic 
errors into account. In addition to the nine QS02s in the 
CDFS, wo also used the emission-line flux data of CXO 
J084837.9-F445352 (CXO 52), a QS02 at z = 3.288 found 
by Stern et al. (2002). Since the flux error for NvA1240 is 
not given by Stern et al. (2002). we assigned a 20% error 
for this line as for the HzRGs mentioned above. Therefore 
the number of QS02s is 10 in total. 

(3) Sy2s: We compiled the flux data of 9 Sy2s observed 

by lUE from the literature. The object name, the compiled 
flux data and their errors, and the reference are given in 
Table 2. For what concerns the data obtained by Thuan 
(1984), note that they reported UV emission- line fluxes 
also for IC 5135, not only for NCG 5135; however, we do 
not include the data of IC 5135 because the signal-to-noise 
ratio of the IC 5135 data is too low (see Thuan 1984 for de- 
tails). For fluxes with no information on their errors in the 
reference, we assigned a 20% error as for the HzRGs and 
QS02s (except for low-accuracy measurements explicitly 
mentioned in the reference for which we assigned a 30% 
error). All of the targets are at z < 0.03 and thus wo re- 
gard this sample as representative of the NLRs in the local 
universe. 



^ http: / /www.mpe. mpg.de/CDFS /data/ 
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Table 2. Emission-Line Measurements of Sy2s 



Object Redshift LyaA1216'' NvA1240=^ CivA1549=^ HeiiA1640=^ Ciii]A1909'' Ref.'' 



NGC 1068 


0.004 


713±100 


224±41 


520±80 


187±29 


240±35 


1 


JNGO 4507 


0.012 


75.6±15.1 


5.2±1.0 


13.5±2.7 


5.6±1.1 


5.8±1.2 


2 


NGC SISS"^ 


0.014 


59.0±11.8 


1.1±0.2 


4.1±0.8 


10.0±2.0 




3 


NGC 5506"^ 


0.006 






4.5±1.4'^ 


2.0±0.6'^ 


3.6±0.7 


2 


NGC 7674 


0.029 


47.0±20.3 




11.4±3.3 


5.1±1.5 


7.9±2.7 


4 


Mrk 3 


0.014 


64±19'* 


Sil"* 


21±2 


9±1 


9±1 


5 


Mrk 573 


0.017 


151.2±22.7 


6.3±0.9 


29.0±4.3 


12.6±1.9 


8.8±1.3 


6 


Mrk 1388 


0.021 






8.3±1.2 


3.8±0.6 


3.6±0.5 


6 


MCG -3-34-64'= 


0.017 


56±11 


5±1 


14±3 


10±2 


7±1 


7 



Fluxes are in units of 10 ergs s ^ cm ^. 
^ References. — (1) Snijders et al. 1986, (2) Bergeron et al. 1981, (3) Thuan 1984, (4) Kraemer et al. 1994, (5) Malkan & Oke 
1983, (6) McAlpine 1988, (7) De Robertis et al. 1988. 
20% of the line fiux is assigned as the flux error. 
30% of the line flux is assigned as the flux error. 



The redshift distributions of the HzRG and the QS02 
samples are shown in Figure 1. The compiled data are 
not corrected for Galactic and intrinsic dust reddening. 
Possible effects of the reddening on our analysis are dis- 
cussed later (§§2.2). Since most of these data were ob- 
tained in low-dispcrsion spectroscopy, the measured Cm] 
flux may be contaminated by the Sini]AA1883,1892 flux. 
However, the contribution of Sim] is thought to be small 
and we will discuss this issue further in §§3.2. 

2.2. Compilation results 

In Figure 2a the compiled line flux ratios are plotted on 
the CivA1549/HeiiA1640 versus Ciii]A1909/CivA1549 di- 
agram. The ClvA1549/HellA1640 flux ratio is expected 
to be sensitive to the gas metallicity. This is because the 
gas temperature decreases when the metallicity increase 
in low-density 10^ cm~^) gas clouds (e.g., Ferland 
et al. 1984) and thus the coUisional excitation of CiV 
is gradually suppressed, while the HellA1640 luminosity 
is basically proportional to the volume of the doubly- 
ionized helium region and thus rather insensitive to the 
gas metalHcity. The Ciii]A1909/CivA1549 ratio is instead 
expected to be sensitive to the ionization degree of gas 
clouds. Therefore, Ciii]A1909/CivA1549 can be used to 
check any dependence of CivA1549/HeiiA1640 on the ion- 
ization state of the gas. Summarizing, a diagnostic dia- 
gram that consists of these two flux ratios is expected to 
be useful to estimate the properties of NLRs only with 
strong UV emission lines, as discussed more extensively 
in §3 and §4 (see also Groves et al. 2004). 

In Figure 2 we also show the effect of the dust extinc- 
tion on the line ratios, for the case of Ay = 5.0 mag. 
An extinction curve described by Cardelli et al. (1989) is 
adopted. Since the dust extinction in typical type 2 AGNs 
is generally Ay ^ 3 mag (see, e.g.. Figure 5 of Nagao et 
al. 2001b), we conclude that the effect of dust extinction 
on our analysis and discussion is not important. 



As clearly shown in Figure 2a, there is no signif- 
icant difference in these two flux ratios between the 
high-z QS02 sample and the local Sy2 sample. Some 
HzRGs show similar flux ratios to QS02s and Sy2s, al- 
though other HzRGs show lower CivA1549/HeiiA1640 and 
higher Ciii]A1909/CivA1549 than QS02s and Sy2s. The 
logarithmically averaged ratios for these three popula- 
tions (excluding upper-limits and lower-limits) are sum- 
marized in Tabic 3. These averaged flux ratios are also 
plotted in Figure 2b. The average flux ratios of HzRGs 
appears systematically different from those of QS02s 
and Sy2s. To sec the statistical significance of the dif- 
ferences in the line flux ratios between HzRGs and 
QS02s - Sy2s, the Kolmogorov-Smirnov (K-S) statisti- 
cal test is applied to the data, discarding upper-limits 
and lower-limits. The null hypothesis is that the flux ra- 
tios (CivA1549/HeiiA1640 and Ciii]A1909/CivA1549) of 
HzRGs (A^data = 49) and QS02s - Sy2s (Adata = 14) 
come from the same underlying population. The derived 
K-S probabilities are 4.0 x IQ-^ for CivA1549/HeiiA1640 
and 8.0 X lO^^ for Cm]A1909/CivA1549. These results 
suggest that the difference in the CivA1549/HeilA1640 ra- 
tio is statistically significant while the difference in the 
Ciii]A1909/CivA1549 ratio is statistically marginal. 

When focusing only on HzRGs, an anti-correlation be- 
tween the emission- line fiux ratios of CivA1549/HeiiA1640 
and Cm]A1909/CivA1549 is apparent in Figure 2a. To 
examine the statistical significance of this correlation, 
the Spearman rank-order test is applied to the HzRG 
data, discarding upper-limits and lower-limits. The de- 
rived Spearman rank-order correlation coefficient (rs) and 
its statistical significance p{rs), which is the probability 
of the data being consistent with the null hypothesis that 
the fiux ratios are not correlated, are rs = —0.54 and 
p{rs) = 2.1 X 10-"^ (Adata = 49). This indicates that the 
apparent anti-correlation between ClvA1549/HellA1640 
and Ciii]A1909/CivA1549 of HzRGs is statistically signif- 
icant. This significance does not disappear by including 
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• HzRG 

O QS02 
A Sy2 




• HzRG 

O QS02 
A Sy2 



Q) 2 



> 
o 
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0.5 - 



Ay = 5 mag 



(b) 



0.2 



0.1 



10 



C III] 



C IV 



Fig. 2. (a) Compiled data plotted on a diagram of 
CivA1549/HeiiA1640 versus NvA1240/CivA1549. Arrows 
denote the upper-limit or lower-limit on the flux ratios. 
Filled circles and arrows denote the HzRG data and open 
circles and arrows denote the QS02 data. Open trian- 
gles denote the Sy2 data. The arrow at the lower-left 
corner in the panel denote the effect of a correction for 
dust extinction {Ay = 5 mag), (b) Logarithmically av- 
eraged flux ratios of HzRGs (filled circle), QS02s (open 
circle) and Sy2s (open triangle) plotted on a diagram of 
CivA1549/HeiiA1640 versus NvA1240/CivA1549. The ar- 
row at the lower-left corner is the same as in (a). Thin 
bars denote the RMS's of the data distribution, and thick 
bars denote the estimated errors on the averaged values. 



QS02s and Sy2s: the corresponding Spearman rank-order 
correlation coefficient and its statistical significance be- 
come rs = -0.47 and p(rs) = 2.3 x 10"^ (/^data = 63). 
This anti-correlation has also boon shown by Allen ct al. 
(1998) for a smaller HzRG sample. We will discuss the 
interpretation of this trend in §§4.1. 



100 



10 



c 1 



0.01 




0.1 



10 100 1000 
Photon Energy (Ryd) 



10-" 



Fig. 3. Adopted SEDs for our photoionization model cal- 
culations. Thin solid line denotes a SED with a strong 
UV bump and thick dashed line denotes a SED with a 
weak UV bump. Thin dotted line denotes the SED used 
by Mathews & Ferland (1983). All the three SEDs are 
normalized to the flux at 1 keV (~73 Ryd). 

Table 3. Averaged Diagnostic Flux Ratios 



Sample 


Number'" 


C IV / He H 


C HI] / C IV 


HzRG 
QS02 
Sy2 


49 
6 
8 


.'-■^^-0.08 


'-'•"'-0.10 
4Q+0.06 



Number of objects for which both of the flux ratios of 
CivA1549/HeiiA1640 and Ciii]A1909/CivA1549 were mesr 
sured. Objects with a lower-limit or upper-limit data are 
not included. 



3. Photoionization models 

To provide a quantitative interpretation of Figure 2, 
we carried out photoionization model calculations. The 
method and the results are given below. 



3.1. Model method 

We performed model runs by using the public photoion- 
ization code Cloudy version 94^ (Ferland 1997; Ferland 
et al. 1998). For simplicity, we assumed uniform gas den- 
sity clouds with a plane-parallel geometry, and we exam- 
ined both dust- free and dusty cases. For the models with 
dust, Orion- type graphite and silicate grains (Baldwin et 
al. 1991; Ferland 1997) were included. Note that the cal- 
culations are not sensitive to the assumption on the grain 
composition (§§4.1.2). Effects of dust grains on the deple- 
tion of heavy elements and on the radiative transfer were 
consistently treated by Cloudy. The parameters for the 
calculations are (1) the spectral energy distribution (SED) 



^ We confirmed for some models that the results of the cal- 
culations do not change significantly if using Cloudy version 
96 instead of version 94: the difference is ~ 10% at maximum. 
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SED with Large UV Bump SED with Small UV Bump 




0.1 0.2 0.5 1 2 5 0.1 0.2 0.5 1 2 

C III] / C IV 



Fig. 4. Model results plotted on a diagram of CivA1549/HeiiA1640 versus Ciii]A1909/CivA1549. Models for gas clouds 
with riH = 10^ cm~^ (thin lines) and 10^ cm~^ (thick lines) are presented. Constant- Zgas sequences and constant-?/ 
sequences are denoted by dotted and solid lines, respectively, (a) Dust-frcc models adopting the input SED with a 
large UV bump, (b) Dusty models adopting the input SED with a large UV bump, (c) Dust-free models adopting the 
input SED with a small UV bump, (d) Dusty models adopting the input SED with a small UV bump. 



of the photoionizing continuum radiation; (2) the hydro- 
gen density of a cloud (nn); (3) the ionization parameter 
{U), i.e., the ratio of the ionizing photon density to the 
hydrogen density at the irradiated surface of a cloud; (4) 
the column density of a cloud (iVn); and (5) the elemental 
composition of the gas. 

As for the SED of the ionizing photons, two extreme 
cases of SED were investigated. The first one is a SED 
with a weak UV thermal bump, which matches the HST 
QSO templates (Zheng et al. 1998; Teller et al. 2002; see 
Marconi ct al. 2004 for more details). The second one 
has a strong UV thermal bump to match the QSO tem- 
plate by Scott et al. (2004). Both SEDs have the same 
optical to X-ray ratio aox (Zamorani et al. 1981), i.e., 
aox = —1.49 (see Elvis et al. 2002), but different slopes 
in the energy range of 9.1eV < hv < 35.5eV; a = —2.17 
for the SED with a strong UV bump and a = —1.71 for 
the SED with a weaker UV bump, where oc v°'. See 
Figure 3 for a graphical representation of the two SEDs. 
Just for the readers' convenience, in Figure 3 we also show 
the SED deduced by Mathews fc Ferland (1983) that has 
been sometimes used for photoionization model calcula- 
tions (see Ferland 1997), although we do not use this SED 
in this work. 



We investigated gas clouds with gas densities njj = 10^ 
cm""^, 10"^ cm~^, 10^ cm~'^ and lO*' cm~^, and ionization 
parameters U = 10~^-^ - 10~°-*, as presented in §§3.2. 
The column density Nn was determined by the criterion 
that the calculations for dust-free gas clouds were stopped 
when the temperature falls to lOOK, below which the gas 
does not c;ontributc significantly to the observed optical 
emission line spectra. Although this lower-limit tempera- 
ture is much lower than other calculations for ionization- 
bounded clouds in the literature, this criterion is neces- 
sary to calculate low-density dust-free gas clouds with a 
high metallicity because the equilibrium temperature of 
such gas clouds is sometimes lower than lOOOK, as it will 
be shown later (see also Ferland et al. 1984). For mod- 
els with dust grains, the stopping temperature was set to 
2000K. This is because the gas tc^iiipc^rature does not de- 
crease efficiently in partially-ionized regions when clouds 
contain dust, due mainly to the depletion of heavy ele- 
ments (i.e., coolants) and to the heating by photoelec- 
trons emitted from grains (see, e.g.. Shields & Kennicutt 
1995 for details). Note that the results of our calculations 
are not sensitive to the lower-limit temperature because 
we are concerned only in relatively high-ionization emis- 
sion lines, which arise in fully- ionized regions, and not 
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within partially-ionized regions. For the chemical com- 
position of dust-free gas clouds, we assumed that the all 
metals scale by keeping solar proportions except for nitro- 
gen, that scales as the square power of other metal abun- 
dances, because nitrogen is a secondary element (see, e.g., 
Hamanii et al. 2002). Here the solar elemental abundances 
are taken from Grevesse & Anders (1989) with extensions 
by Grevesse & Noels (1993). For dusty gas clouds, we as- 
sumed the depleted gas-phase elemental abundance ratios 
by adopting the depletion factors given by Ferland (1997). 

3.2. Model results 

In Figures 4, 5 and 6, the results of the model calculations 
are plotted on a diagram with CivA1549/HeiiA1640 and 
Ciii]A1909/CivA1549, for both dust-free and dusty cases. 
Models for low-density gas clouds (nn = 10^ cm~^ and 
10^ cm~^) are presented in Figure 4, those for high-density 
gas clouds (nn = 10'' cm^^) are presented in Figure 5, and 
those for very high-density gas clouds (nn = 10^ cm~^) 
are presented in Figure 6. For low-density cases, the dif- 
ference of the model results between those with nu = 10^ 
cm~^ and those with nn = 10^ cm~^ is very small. 

For low-density gas clouds without dust grains, the 
flux ratio of CivA1549/HeiiA1640 strongly depends on the 
gas metalUcity while Ciii]A1909/CivA1549 allows to re- 



move the degeneracy with U, and thus the diagram with 
CivA1549/HeiiA1640 and Ciii]A1909/CivA1549 is a good 
metallicity diagnostics (Figures 4a and 4c). The strong 
variation of CivA1549/HeiiA1640 with metallicity is due 
to the strong dependence of ClvA1549 emissivity on the 
gas temperature, which decreases rapidly with metallicity 
(in dust-free clouds) due to an efficient cooling by infrared 
fine-structure lines (e.g., Ferland et al. 1984). The metal- 
licity dependences of some infrared fine-structure lines 
([Oiii]88/um, [Niii]57/im and [Neiii]15.6/um) and gas tem- 
perature for gas clouds with nu = 10^ cm~^ are shown in 
Figure 7. Here we adopt the averaged temperature over 
doubly-ionized helium regions as a representative temper- 
ature in photoionization equilibrium. The fluxes of the 
fine-structure lines plotted in Figure 7 are normalized to 
the H/3 flux, because the H/3 flux scales with the ionization 
photon flux. 

The metallicity dependence of CivA1549/HeiiA1640 
in dust-free gas clouds becomes smaller at higher gas 
densities (Figures 5a, 5c, 6a, and 6c). This is mainly 
due to the decreased cooling efficiency of infrared fine- 
structure lines, which are suppressed in high-density 
gas clouds owing to collisional de-excitation (Figure 7). 
However, even when gas clouds with nn — 10^ cm""^ 
are concerned, the diagram with CivA1549/HeiiA1640 and 
Ciii]A1909/CivA1549 maybe useful to constrain the gas 
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metallicity, because the ClvA1549/HeiiA1640 ratio de- 
pends on metallicity even at nn = 10^ cm~^, although the 
inferred metallicity is different depending on the specific 
gas density (Figures 5a and 5c). The CivA1549/HellA1640 
versus Ciii]A1909/CivA1549 diagram is no more useful 
when the gas density is very high, njj = 10^ cm^'^ 
(Figures 6a and 6c). The metallicity dependence of the 
CivA1549/HciiA1640 has almost disappeared at this gas 
density. However, as shown later, the very high-density 
models (uh = 10® cm~^) do not provide a good descrip- 
tion of the observed data. 

Gas clouds with dust grains show only a very small de- 
pendence of CivA1549/HeiiA1640 on metallicity (Figures 
4b, 4d, 5b, 5d, 6b, and 6d). This is mainly because 
the equilibrium temperature of gas clouds does not drop 
off significantly when gas metallicity is high, thanks to 
the presence of dust grains (Figure 7d; sec Shields & 
Kennicutt 1995 for more details). This result is almost 
independent of the adopted SED, gas density, and ioniza- 
tion parameter. 

Figures 4, 5 and 6 indicate that the results with a 
large UV bump SED and with a small UV bump SED are 
similar. We thus conclude that SED effects on our analysis 
and discussion are negligible. In the following discussion, 
we only refer to the models with a small UV bump SED. 

Note that the diagnostic diagram on which we are 
focusing has been investigated for various purposes in 
the past, since all the three lines (HeiiA1640, CivA1549 



and Clil]A1909) are very strong and the wavelength sep- 
aration is small (i.e., their ratios arc little sensitive to 
dust extinction). Allen ct al. (1998) investigated pho- 
toionization models and fast-shock models showing that 
the data of HzRGs on the CivA1549/HeiiA1640 versus 
Ciii]A1909/CivA1549 diagram are consistent with pho- 
toionization, but are hard to be accounted for by fast- 
shock models. Groves et al. (2004) showed that this dia- 
gram is useful to examine the gas metallicity of NLRs and 
mentioned that HzRGs may have gas clouds with sub-solar 
metallicity. We use this diagram to analyze the gas metal- 
licity of NLRs for a larger sample of high- 2; narrow-line 
AGNs and to investigate the possible metallicity evolu- 
tion of NLRs. 

Finally we tackle the issue of whether the flux of 
Ciii]A1909 is contaminated by the Siiii]AA1883,1892 
emission significantly. The predicted ratio of 
Siiii]AA1883,1892/Ciii]A1909 is plotted as a function 
of ionization parameter in Figure 8. Here we investigate 
gas clouds with (riR, ^gas) = (10^ cm'^, 0.5 Zq), (10^ 
cm-3, 2.0 Zq), (IQS cm-3, 0.5 Zq), and (10^ cm'^, 2.0 
Zq). For clouds with dust grains, the predicted ratio is 
~0.01 regardless of density, metallicity, and ionization 
parameter. The contribution of Silll]AA1883,1892 is 
thus negligible when dusty clouds are concerned. For 
dust-frcc clouds, the predicted ratios arc higher, i.e., 
~0.1. However, this is still significantly smaller than the 
measurement errors on Cili]A1909 fluxes. Therefore, we 
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conclude that the contribution of Siiii]AA1883,1892 does 
not affect our results and discussion signiflcantly. 

4. Discussion 

4.1. Comparison of the data with models 
4.1.1. Shock models 

Before comparing the data presented in §2 with the pho- 
toionization models presented in §3, we examine whether 
shock models can explain the observed flux ratios. In 
Figure 9, we compare the data of the HzRG, QS02 and 
Sy2 samples with the shock models presented by Dopita 
& Sutherland (1996). Both pure shock-only models and 
shock plus precursor models are examined here: shock- 
only models with a shock velocity of 150 km s~^ < i^shock 

< 500 km s~^ and with a magnetic parameter of /xG 
cm~^/^ < B jn}-!'^ < 4 /xG cm""^/^, and shock plus precur- 
sor models with a shock velocity of 200 km s^^ < Wghock 

< 500 km s~^. As shown in Figure 9, both pure shock- 
only models and shock plus precursor models underpre- 
dict the flux ratio of Ciii]A1909/CivA1549 and thus can- 
not explain the observed data. This suggests that the main 



ionization mechanism of NLRs (or at least the CivA1549, 
HellA1640 and Ciii]A1909 emitting regions in NLRs) is not 
associated with shocks, but is photoionization. Figure 9 
also suggests that the difference in the flux ratios between 
HzRGs and QS02s Sy2s cannot be ascribed to shocks. 
Instead, a preferential contribution of shocks to QS02s 
(not to HzRGs) is required if shocks are at the origin of 
the difference in the flux ratios. This requirement is in 
the opposite direction of the natural expectation (HzRGs 
should be more affected by shocks because of jets and ex- 
panding radio lobes). Therefore, the systematic difference 
in the flux ratios between HzRGs and QS02s - Sy2s must 
be ascribed to causes other than shocks. This issue will be 
discussed further in §§4.1.3. 

4.1.2. Dust grains 

In Figure 10, we compare the observational data with the 
results of our photoionization model calculations, both 
with and without dust grains. Contrary to the dust-free 
models, the dusty models predict too narrow ranges of the 
CivA1549/HeiiA1640 flux ratio, regardless of the gas den- 
sity (Figures 10b, lOd, and lOf). The CivA1549/HeiiA1640 
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and Zgas = 2.0Zq, respectively, (a) Dust-free models with 
a large UV bump SED. (b) Dusty models with a large UV 
bump SED. (c) Dust-frce models with a small UV bump 
SED. (d) Dusty models with a small UV bump SED. 



flux ratio varies only in a factor of 3 at maximum even 
when the metallicity varies in the range 0.2 < Z^^s/Zq < 
10.0. More importantly, the dusty models cannot explain 
values of CivA1549/HeiiA1640 larger than ~ 1.5, which 
are instead observed in most sources. These results sug- 
gest that the CivA1549, HeiiA1640 and Ciii]A1909 emit- 
ting regions in NLRs are dust-free. This is consistent with 
previous works that gas clouds in the high-ionization part 
of NLRs are dust-free (e.g., Marconi et al. 1994; Ferguson 
et al. 1997; Nagao et al. 2003). We verified that the effects 
of changing the grain composition (see Ferland 1997 for 
details) are less than 30% on the line flux ratios. We thus 
conclude that the models with dust grains are not appro- 
priate to interpret the ratios of ClvA1549/HellA1640 and 
Ciii]A1909/CivA1549. 



4.1.3. Ionization parameter 
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Ciii]A1909/CivA1549 (Figures 10a and 10c), while the 
very high-density models (nn = 10® cm~^) cannot ex- 
plain most of the observational data (Figure lOe). And we 
can deduce that uh = 10^ cm""^ is an upper bound to 
the average NLR density. Both high-density models and 
low-density models suggest that the observational data 
are consistent with clouds with an ionization parame- 
ter 10-2-2 ^ U ^ 10-1-^ for HzRGs and U w lO'^-^ 
for QS02s and Sy2s (Figures 10a and 10c). The model 
sequences with a constant ionization parameter have 
negative slope in the diagram of CivA1549/HeiiA1640 
versus Ciii]A1909/CivA1549, which is consistent with 
the observational trend described in §§2.2. Therefore 
the anti-correlation between CivA1549/HeiiA1640 and 
Ciii]A1909/CivA1549 seen in HzRGs can be naturally 
explained by the wide range of the ionization parame- 
ter for HzRGs. Although the inferred range of ionization 
parameter is higher than the range estimated through 
rest-frame optical diagnostics (which use lower ionization 
emission lines) reported by some previous works (e.g., 
Ferland & Netzer 1983; Ho et al. 1993), it is consistent 
with previous estimates of the ionization parameter for 
the high-ionization parts in NLRs (e.g., Oliva et al. 1994; 
Murayama & Taniguchi 1998; Nagao et al. 2001a, 2001b). 



4.1.4. Gas density and metallicity 



density (nn = 10 cm '') can successfully explain the As mentioned in the last subsection, photoionization 



observed range of flux ratios CivA1549/HeiiA1640 and models with nn 



10° cm predict too narrow 
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ranges of the flux ratios of CivA1549/HeiiA1640 and 

Ciii]A1909/CivA1549 with respect to the observed ranges. 
This suggests that the typical density of the CivA1549, 
HellA1640 and Clll]A1909-cniitting region in the NLR 
should be less than 10^' cm^'\ Therefore, in the follow- 
ing discussion, we consider only the two models, i.e., the 
low-density models with nn = 10^ cm~^ (note that the 
results are very similar if adopting nn = 10^ cm~^) and 
the high-density models with riu = 10^ cm~^. 

As shown in Figures 10a and 10c, the estimated metal- 
licity is different when different gas densities are adopted. 
Models with nn = 10^ cm"'' suggest sub-solar metal- 
licities (0.2 ^ Zge.s/ZQ 1.0; Figure 10a) while mod- 
els with nu = 10^ cm""' siiggest a wide metallicity 
range (0.2 ^ Zg^/ZQ ^ 5.0; Figure 10c). The "low- 
density and sub-solar metallicity" scenario appears to be 
consistent with the results reported by Iwamuro et al. 
(2003), while the "high-density with wide metallicity dis- 
tribution" scenario appears in agreement with the re- 
sults reported by De Breuck et al. (2000). Although some 
observational data with lower CivA1549/HeiiA1640 and 
Ciii]A1909/CivA1549 ratios appear to deviate from the 
model predictions with nu = 10^ cm"'' (Figure 10c), 
these deviations can be reconciled by introducing a mod- 
erate amount of extinction. Although the above two sce- 



narios are hard to be discriminated only by means of 

the CivA1549/HcnA1640 versus Ciii]A1909/CivA1549 di- 
agram, the actual situation may be intermediate between 
the two possible scenarios (§§4.2). Note that the abso- 
lute values of the inferred gas metallicity are not very 
accurate due to a density dependence of the flux ratio 
CivA1549/HeiiA1640 at high-densities (nu 10^ cm-^). 
Nevertheless the diagnostic diagram in Figure 10 is use- 
ful if we are interested in the relative metallicity trends 
of NLR clouds, or in the evolution of gas metallicity, 
adopting the assumption that there are no strong corre- 
lations between gas density and redshift. Finally, we note 
that the density affects the inferred metallicity mostly at 
Zgas 0.5^0. At low mctallicities (^gas ^ 0.2Zq) the 
models are less sensitive to the gas density (Figures 10a 
and 10c). This result is important, since it allows to use 
this diagnostic diagram to identify low metallicity NLRs, 
regardless of the gas density. 

4.2. Comparison with the NvX1240 diagnostics 

The gas metallicity of the BLR clouds in QSOs has 
been often investigated by using the diagnostic flux ra- 
tios of NvA1240/CivA1549 and NvA1240/HeiiA1640 (e.g., 
Hamann & Ferland 1992, 1999; Dietrich et al. 2003; 
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versus CivA1549/HeiiA1640 (d). 



Nagao ct al. 2005). This method has been sometimes 
applied also to the NLR clouds (e.g., van Ojik et al. 
1994; Villar-Martin et al. 1999; De Breuck et al. 2000; 
Vcrnct ct al. 2001). Therefore it is interesting to com- 
pare our diagnostics with the NvA1240 diagnostics. In 
Figure 11, the observational data are compared with 
the photoionization models for nn ~ 10^ cm~^ and 
riH = 10^ cm-3 on the diagram of NvA1240/CivA1549 
versus CivA1549/HeiiA1640, and NvA1240/HeiiA1640 
versus CivA1549/HciiA1640. Although high ratios of 
NvA1240/CivA1549 and NvA1240/HeiiA1640 are some- 
times interpreted simply as a result of high metallici- 
tics, the low-density models (tt-h = 10^ cm~^) with a 
high metallicity predict not only high NvA1240/CivA1549 
ratios but also very low CivA1549/HeiiA1640 ratios, 
in disagreement with the observed data (Figure 11a). 
This is mainly attributed to the decreased equilib- 
rium gas temperature due to the efficient cooling by 
infrared fine-structure lines as investigated in §§3.2 
(Figure 7). Moreover the low-density models predict 
NvA1240/HeiiA1640 < 0.4 when U < 10"°-^ is concerned, 
which is also in disagreement with the observed data 
(Figure lib). The high-density models (tzh = 10^ cm~^), 
on the other hand, predict higher CivA1549/HeiiA1640 



ratios than the low-density models owing to the sup- 
pressed cooling. Since the high-density models predict 
higher CivA1549/HeiiA1640 ratios than the observed val- 
ues (Figures 11c and lid), the typical gas density of (the 
high-ionization part of) the NLRs may be lower than 
riH = 10^ cm~^. 

Although the NvA1240-detected objects appear to 
have super-solar metallicity, the NvA1240 fluxes of more 
than half of the objects in the sample are upper-limits 
and they might be consistent with sub-solar metallicities. 
Summarizing, since NvA1240 becomes very weak (gener- 
ally undetected) for Zg^s ^ '2Zq NLRs, it is very difficult 
to investigate the metallicity evolution of the NLR based 
on the NvA1240 diagnostics. 

4.3. Metallicity evolution of the NLRs in AGNs 

Now we discuss whether the gas metallicity of the NLRs 
in high-z AGNs evolves as a function of rcdshift. We first 
focus only on the HzRG data, because HzRGs and QS02s 
may have a different NLR properties, as discussed above. 
To investigate this issue, we divide the sample of ob- 
jects having a measurement of both CivA1549/HeiiA1640 
and Ciii]A1909/CivA1549 (49 objects, discarding the ob- 
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jects with upper-limits or lower-limits) into three groups: 

1.2 < z < 2.0 (17 objects), 2.0 < z < 2.5 (20 objects), 
and 2.5 < z < 3.8 (12 objects). Note that most of the 
highest- 2; sample are at 2.5 < z < 3.0 and only two ob- 
jects are at z > 3 (Figure la). Wc calculate the logarith- 
mically averaged flux ratios of CivA1549/HeiiA1640 and 
Ciii]A1909/CivA1549 for these three groups. The results 
arc given in Tabic 4 and plotted on the diagnostic dia- 
gram in Figures 12a and 12c. The thick errorbars give the 
estimated errors on the means. As clearly seen, no evolu- 
tionary tendency in the flux ratios within our HzRG sam- 
ple is found. The higher- 2; objects do not show significant 
metallicity decrease with respect to the lower- 2; objects, 
at variance with the result reported by Dc Breuck et al. 
(2000). Although the absolute value of the gas metallicity 
derived from Figure 10 is subject to non-negligible un- 
certainties, as discussed in §§4.1.4, our conclusion is not 
affected by this issue in terms of relative gas metallicities, 
i.e., the gas metallicity in NLRs of HzRGs does not change 
significantly in the rcdshift range 1.2 < z < 3.8, or more 
conservatively, in the range 1.2 < 2 < 3 (owing to the lack 
of objects at 2; > 3). 

Recently Nagao et al. (2005) reported that the gas 
metallicity of the BLR in QSOs with a given luminos- 
ity is independent of rcdshift in the range 2.0 < z < 4.5. 
The latter result is consistent with that obtained for the 
NLR metallicity of HzRGs presented in this paper. Nagao 
et al. (2005) also reported that the BLR metallicity is 
tightly correlated with the QSO luminosity (see also, e.g., 
Hamann & Ferland 1993, 1999). Motivated by this corre- 
lation between BLR metallicity and QSO luminosity, we 
have investigated whether the NLR metallicity is corre- 
lated with the AGN luminosity or not. However, it is very 
difficult to measure the luminosity of both AGNs and their 
host galaxies when type-2 AGNs are concerned. This is 
because the central engine is hidden by the dusty torus 
and because the broad-band photometric flux is largely 
attributed to the nebular emission, not only to the stel- 
lar continuum emission. We thus adopt the HciiA1640 
emission-line luminosity [i(Hell)], as an indicator of the 
AGN luminosity. This assumption is based on the fact that 



the HellA1640 luminosity is simply proportional to the 

volume of the doubly-ionized helium region, which scales 
to the AGN luminosity. We calculate L(Heii) from the 
HeiiA1640 flux by adopting a cosmology with (Otot, f^m, 
Oa) = (1.0, 0.3, 0.7) and Hq = 70 km s'^ cm"! Mpc"!. 
L(Heii) is not corrected for the slit loss, which may be 
non-negligible for some cases. To investigate the depen- 
dence of the line flux ratios on L(Heii), we divide HzRGs 
into three groups: 41.5 < log L(Heii) < 42.5 (13 objects), 
42.5 < log L(Heii) < 43.0 (21 objects), and 43.0 < log 
L(Heii) < 45.0 (15 objects), where L(Heii) is in imits of 
erg s~^. The logarithmically averaged flux ratios and the 
RMS's of CivA1549/HeiiA1640 and Ciii]A1909/CivA1549 
for these three groups are given in Table 4, and are plot- 
ted on the diagnostic diagram in Figure 12b and 12d. 
The thick errorbars give the errors on the means. As 
a result, we find a systematic trend in our HzRG sam- 
ple that the HzRGs with larger L(Heii) (i.e., more lu- 
minous AGNs) tend to have lower CivA1549/HeiiA1640 
and Cili]A1909/CivA1549 flux ratios. This result is con- 
sistent to the interpretation that the NLR in more lu- 
minous HzRGs have higher metallicity gas clouds. This 
"luminosity-metallicity relation" for the NLR in HzRGs 
is in agreement with the same relation seen for the BLRs 
in high-2 QSOs. 

The similarity of the flux ratios between the QS02 
sample and the Sy2 sample may suggest that the gas 
metallicity of non-radio-selected AGNs does not evolve 
signiflcantly from z ~ 3 to the local universe. However 
it should be kept in mind that there is a large differ- 
ence in the absolute luminosity between the Sy2 sam- 
ple and the QS02 sample. This result should be inter- 
preted in the sense that the NLR metallicity is not sig- 
niflcantly different between faint local AGNs and bright 
high-2 (1.5 ^ z ^ 3.7) AGNs. If the luminosity-metallicity 
relation discussed above is taken into account, this result 
might imply a metallicity evolution of the NLR with a 
given luminosity from high-z to z ~ 0. Surveys and spec- 
troscopic studies on type-2 faint AGNs at high-z are re- 
quired to investigate this issue further. 
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Fig. 12. Averaged flux ratios compared with the model predictions on the CivA1549/HeiiA1640 versus 
Ciii]A1909/CivA1549 diagram. Filled circle, open circle and open triangle denote the averaged values of the HzRG 
sample, the QS02 sample and Sy2s sample. Models with nn = 10^ cm~^ arc plotted in panels (a) and (b), while 
models with uh = 10^ cm~"^ are plotted in panels (c) and (d). For both cases, models with a small UV bump SED and 
without dust are plotted. The redshift dependence of the flux ratios of the HzRGs are investigated in panels (a) and 
(c), where filled diamond, square and star denote the averaged flux ratios of the HzRGs at 1.2 < z < 2.0, 2.0 < z < 2.5 
and 2.5 < z < 3.8, respectively. The luminosity dependence of the flux ratios of the HzRGs are investigated in panels 
(b) and (d), where filled diamond, square and star denote the averaged flux ratios of the HzRGs with 41.5 < log 
L(Heii) < 42.5, 42.5 < log L(Heii) < 43.0, and 43.0 < log L(Heii) < 45.0 respectively, where L(Heii) is in units of 
ergs s~^. Error bars are as in Figure 2b. 



It has been observationally confirmed that QSOs in- 
cluding ratio galaxies tend to reside in massive elliptical 
galaxies at least in the local universe (e.g., McLure et al. 
1999; Dunlop et al. 2003; Floyd et al. 2004) and proba- 
bly also in the high-z universe (e.g., Falomo et al. 2005; 
Kuhlbrodt et al. 2005). Our results suggest that the main 
star-formation event in such massive host galaxies is com- 
pleted. This picture is consistent with the recent findings 
of massive evolved galaxies at z > 1.5 (e.g., Cimatti et al. 
2004; McCarthy et al. 2004; Labbe et al. 2005; Saracco et 
al. 2005). The non-evolution of the gas metallicity of the 
NLR in our sample suggests that the epoch of major chem- 
ical enrichment in the host galaxies of QS02s and HzRGs 
must have occurred at z ^ 3. If the minimum timescalc 
for significant enrichment of carbon (~0.5 Gyr) is taken 
into account, the major epoch of star formation may be 
at .2 ^ 4. 

Finally we briefly discuss a specific, interesting QS02, 
CDFS-901 (Figure 2a). Although the emission-line spec- 



trum of this object is hard to be explained by photoion- 
ization models as seen in Figure 10, it appears to be con- 
sistent with pure shock-only models as shown in Figure 
9. This may suggest that the NLR of CDFS-901 is dom- 
inated by shocks. However, there is another possibility, 
that is CDFS-901 may be a narrow-line type-1 QSO, 
i.e., a brighter analogue of narrow-line Seyfert 1 galax- 
ies (NLSls). If the emission-line spectrum of CDFS-901 
is not from its NLR but from its BLR, its very large 
CivA1549/HeiiA1640 ratio (> 7.0) is naturally explained, 
since this flux ratio is expected to be 10 for BLRs (al- 
though its Ciii]A1909/CivA1549 ratio ~0.2 is very smafl 
for a BLR, e.g., Nagao et al. 2005). Since the X-ray spec- 
tral slope (hardness ratio) of this object was not measured 
due to the lack of photon statistics, its nature remains am- 
biguous based on the currently available data. This object 
is interesting because NLSls are sometimes thought to be 
AGNs with super-massive black holes in their growing-up 
phase (e.g., Mathur 2000). If this object is a really high- 2; 
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analogue of NLSls, it may be a very interesting target to 
investigate the evolution of AGN activities and supermas- 
sive black holes. 

5. Summary 

In order to investigate the possible metallicity evolution 
of NLR gas clouds in AGNs, we compiled the fluxes of 
CivA1549, HeiiA1640 and Ciii]A1909 for a large sample of 
narrow-line AGNs, including HzRGs, high-z QS02s, and 
local Sy2s. Since all of these three emission lines are mod- 
erately strong even in the faint HzRGs and QS02s, this 
approach enables us to investigate a large number of such 
objects. By comparing the compiled flux ratios with the re- 
sults of our photoionization model calculations, we found 
the following results. 

— The observational data are inconsistent with the pre- 
dictions of shock models, suggesting that the NLRs are 

mainly photoionizcd. 

— The photoionization models with dust grains predict 
too narrow ranges of flux ratios and in disagreement 
with the observed ranges, suggesting that the high- 
ionization part of NLRs (on which we focused in this 
work) is dust-free. 

— The ionization parameter of (the high-ionization parts 
of) NLRs is estimated to be IQ-^ ^ ^ [/ ^ IQ-i '' for 
HzRGs and U « IQ-i ^ for QS02s and Sy2s. 

— The photoionization models with nn = 10^ cm^"* can- 
not explain the observational data, suggesting that the 
typical gas density is lower than 10^ cm~'^. 

— There are two possible interpretations for the observed 
data: low-density gas clouds (nn ^ 10'^ cm^'') with a 
sub-solar metallicity (0.2 ^ Zg.^s/ZQ ^ 1.0), or high- 
density gas clouds (nn ^ lO"' cm^'*) with a wide range 
of gas metallicity (0.2 <, Z^^s/Zq <, 5.0). 

— Our method using only the flux ratios of 
CivA1549/HeiiA1640 and Ciii]A1909/CivA1549 is 
particularly useful to examine relative difference 
in gas metallicity of NLR clouds, or to investigate 
possible metallicity evolution of NLRs, although 
the inferred absolute values of metallicity contain 
non-negligible uncertainties. 

— Regardless of the density, the proposed diagnostic dia- 
gram is useful to identify low metallicity NLRs {Zg^s ~ 
O.2Z0). 

— We find no evidence suggesting a significant evolution 
of the gas metallicity in the NLRs of HzRGs in the 
redshift range 1.2 < z < 3. 

— We find a systematic trend for more luminous AGNs 
to have more metal-rich NLRs (luminosity-metallicity 
relation), which is in agreement with the results from 
the studies on the BLRs. 

— The non-evolution of the gas metallicity of the NLRs 
implies that the major epoch of star formation in the 
host galaxies is at z ^ 4. 
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